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Abstract 


A  performance  analysis  of  the  thermochemical/ 
thermomechanical  behavior  of  graphite  is  made 
extending  present  analysis  procedures  to  include 
the  allotropic  features  of  bulk  graphite.  The  anal¬ 
ysis  comprises  a  complete  definition  of  graphite 
sublimation,  and  formulates  a  model  for  mechan¬ 
ical  erosion  based  on  erosion  of  individual  grains. 
Newly  computed  thermodynamic  functions  for  the 
polyatomic  ring  and  chain  carbon-molecules  are 
used  to  define  the  sublimation  regime;  the  probable 
uncertainty  in  these  functions  on  computed  vapor 
pressure  and  mass  transfer  is  evaluated.  Removal 
of  graphite  filler  grains  owing  to  preferential 
erosion  of  the  binder  matrix  is  described;  numer¬ 
ical  results  are  presented  for  the  fraction  of  solid 
material  mechanically  eroded  from  the  surface. 
Models  fabricated  from  graphites  encompassing 
fine- to  coarse-grain  filler  particles  were  tested 
over  a  wide  environmental  spectrum.  The  high- 
pressure  test  data  are  correlated  to  establish  the 
creditability  of  the  analysis  procedure.)  The  com¬ 
posite  effect  of  thermochemical  and  th^rmomech- 
anical  erosion  is  depicted  in  the  test  data 
correlations. 
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Blowing  value  of  Stanton  number 
Nonblowing  value  of  Stanton  number 
Binary  diffusion  coefficient 
Activation  energy 

Internal  energy  per  interior  carbon  atom 

Weight  fraction 

Standard  free  energy 

Standard  free  energy  of  formation 

Height  of  exposed  filler  particle  prior  to 
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Characteristic  dimension  of  filler  particle 
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referenced  to  absolute  zero  temperature 
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orbitals 
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Standard  heat  of  formation  at  0°K 
Sublimation  energy  of  graphite 

Parameter  in  Eqs.  7  to  10  equal  "j/n; 
j=l,2,...  n-lin  Eq.  7;j  =  2,  3,  ...  n-1  in 
Eq.  8;  and  j=l,  2,  .  .  .  n  in  Eqs.  9  and  10 

Conductivity 

Mass  fraction  of  i^'^  chemical  element 
Mass  fraction  of  j  "  chemical  specie 

Equilibrium  constant  at  constat  pressure 
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Lewis -Semenov  number  -  - 1- 
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Molecular  weight,  Mach  number 
Mass  of  carbon  atom 
Mass  flux 

Number  of  carbon  atoms 


P  Pressure 

q  Heat  transfer  rate 

r  Radial  distance  in  ring  molecule 

R  Universal  gas  constant 

Rq^  Resonance  energy  of  species 

B  Carbon-carbon  distance 

s  Surface  recession  rate 

S  Distance  along  body  surface 

S..  Overlap  integral  of  i^^  and  atomic 

orbitals 

T  Temperature 

u,  V  Components  of  velocity  parallel  and 

normal  to  surface 

X  Moving  coordinate 

y  Fixed  coordinate,  distance  measured  from 

initial  surface  position,  or  boundary  layer 
coordinate 


Greek  Symbols 


c  Emissivity,  energy  of  molecular  orbital 

0  Time 

P  Density 

(7  Stefan- Boltzmann  constant 

Xg  Bending  force  constant 

Xg  Stretching  force  constant 

X.J,  Torsional  force  constant 

bi  Fundamental  vibrational  frequency 
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Subscripts 

B  Bincic  r,  or  bulk 

C  Carbon 

D  Diffusion  controlled  oxidation 

e  Boundary-laye i'  edge 

F  Filler 

Gasi's 

i  Initial,  or  i^^  element 

•  th 

j  j  specie 

L,  Linear 

R  Ring 

s  Solid  .naterial,  or  sublimation 

S  Stagnation 

w  wall 

I.  Introduction 

Carbon  and  its  allotropic  forms  have  been 
used  extensively  for  thermal  protection  of  super¬ 
orbital  and  strategic  entry  vehicles  over  the  past 
decade.  These  materials  are  thermally  attractive 
because  of  their  ability  to  support  high  surface 
temperatures  with  comparatively  low  surface 
erosion.  Although  the  carbonized  and  graphitic- 
carbon  allotropes  have  been  studied  analytically 
with  some  vigor,  recently  developed  high-pressure 
test  facilities  have  revealed  a  large  deficiency  in 
analytical  predictions  of  the  observed  mass-loss 
data.  As  a  consequence,  resort  to  empirical  and 
semi-empirical  correlations  of  the  test  data  has 
ri'sulted  with  rather  drastic  results  evidenced  in 
extrapolation  to  actual  flight  environments. 

Understandably,  empirical  correlations  oftest 
data  are  at  best  provisional;  little  insight  into  the 
actual  mechanisms  operative  during  the  ablation 
process*  s  is  gained^L  2).  Conversely,  rather 
sophisticated  analyses  of  the  thermochemical 
erosion  of  a  variety  of  carbon-based  materials  has 
evolved^^'*^’ 5),  Development  of  these  analyses  has 
been  closely  followed  by  a  number  of  papers  deal¬ 
ing  with  various  aspects  of  the  thermomechanical 
ablation  problem)^* 8),  and  are  specialized  to  the 
particular  material  under  study.  The  composite  of 
these  analyses  has  failed  to  provide  repetitive  cor¬ 
relation  of  the  test  data- -obviating  their  generality. 
In  retrospect,  it  appears  that  the  physical  model 
used  to  describe  the  ablation  mechanisms  has  been 
neglected  in  deference  to  mathematical  eloquence. 

It  is  now  evident  that  a  large  number  of  disciplines 
are  required  to  effectively  approach  the  analytical 
description  of  even  what  previously  was  described 
as  homogeneous  materials;  e.  g.  ,  graphite. 

This  paper  is  directly  concerned  with  the 
ablation  performance  of  graphite  in  high-pressure 
environments.  Here,  the  basic  structure  of  bulk 
graphite  is  examined  and  used  to  formulate  a 
model  to  describe  the  surface  erosion  mechanism. 
In  the  analysis  development,  attention  is  directed 
toward  the  heterogeneous  structure  of  graphite,  the 
formation  of  large  polyatomic  carbon-vapor 
species  subliming  from  the  heated  surface,  and  the 
preferential  erosion  of  dissimilar  constituents 
within  the  graphite  structure  promoting  removal  of 
particulate  solid.  To  achieve  the  analytical  per¬ 
formance  description,  newly  computed  thermo¬ 
dynamic  functions  for  the  polyatomic  ring  and  chain 


carbon  molecules  are  used  to  define  the  sublima¬ 
tion  regime;  the  probable  uncertainty  in  these 
functions  on  computed  vapor  pressure  and  mass 
transfer  is  evaluated.  This  uncertainty  is  carried 
throughout  the  paper  to  assess  the  sensitivity  on 
computed  surface  recession  and  energy  transfer  in 
a  self-consistent  manner.  Removal  of  graphite 
filler  grains  caused  by  preferential  erosion  of  the 
binder  matrix  is  formulated;  numetical  results  are 
presented  for  the  fraction  of  solid  material  eroded 
from  the  surface.  To  establish  the  solid  fraction 
of  eroded  material,  ablation  test  data  in  the  10-  to 
100-atm  pressure  range,  encompassing  fine-  to 
coarse-grain  filler  particles,  are  correlated.  The 
composite  effect  of  thermochemical  and  thermo¬ 
mechanical  erosion  is  depicted  in  the  test  data  cor¬ 
relations.  A  priori,  it  should  be  understood  that 
the  proposed  model  is  provisional  in  that  the  high- 
pressure  test  data  lack  the  resolution  necessaryfor 
detailed  qualification  of  the  analysis  procedures 
(i.  e. ,  the  individual  grains  are  too  small  to  be 
observed  leaving  the  surface). 

II.  Bulk  Graphite 

There  is  an  understandable  tendency  to  treat 
graphite  as  a  homogeneous  solid  when  analyzing 
ablation  processes;  graphite  consists  of  a  single 
chemical  element- -carbon.  Bulk  graphite,  by  the 
nature  of  the  process  used  in  its  formation,  is  far 
from  homogeneous  and  failure  to  account  for  the 
heterogeneous  structure  leads  to  an  incomplete 
description  of  the  surface  erosion  mechanism.  In 
this  paper,  bulk  graphite  denotes  the  different 
forms  of  carbon  existing  in  the  solid  caused  by 
differences  in  constituent  materials  forming  the 
filler  grains  and  binder  matrix. 


Physical  Characteristics 

An  understanding  of  the  microstructure  of 
graphite  requires  some  knowledge  of  the  processes 
by  which  graphite  is  formed.  The  commercial 
graphites  used  in  aerospace  applications  can  be 
described  as  synthetic  or  manufactured.  These 
bulk  graphites  are  manufactured  by  mixing  an 
organic  binder,  usually  a  coal  tar  pitch,  with 
carbon  filler  particles  and  molding  or  extruding  the 
mixture  to  form  a  "green"  billet.  The  binder  is 
then  carbonized  by  heating  to  temperatures  on  the 
order  of  1500'F,  followed  by  further  heating  to 
temperatures  up  to  5400° F,  to  produce  a  more 
ordered  microstructure.  The  filler  material  may 
be  any  particulate  graphite  or  carbon,  such  as 
natural  graphite,  carbon  black,  coke,  etc.  For 
economic  reasons,  filler  material  is  most  com¬ 
monly  made  from  the  coked  residue  of  the  petro¬ 
leum  refining  process.  The  coke  mass  is  crushed 
and  ground  to  eliminate  gross  porosity,  and 
screened  into  controlled  particle  sizes.  Because  of 
unavoidable  intercrystalline  micropores,  the 
particles  retain  a  minimum  porosity  of  5  to 
resulting  in  apparent  densities  of  about  2.  0  to 
2.  2  gm/cc. 

Figure  1  provides  a  schematic  of  the  proposed 
bulk  graphite  arrangement  after  graphitization’t. 

The  filler  particles  appear  as  an  array  of  loosely 
packed  irregular-shaped  coke  elements  cemented 
together  by  the  carbonized  binder.  The  ratio  of 


<'The  description  given  here  abets  the  excellent 
discussion  of  graphite  by  Riley^®®^. 


binder-to-filler  weight  may  be  varied  over  wide 
limita,  with  normal  practice  being  between  20  to 
40%.  In  the  baking  cycle  (1500*F)  the  pitch  binder 
is  pyrolyzed,  losing  about  40  to  50%  of  its  initial 
weight.  During  the  baking  cycle,  and  upon  subse¬ 
quent  graphitization,  the  binder  shrinks.  Since 
shrinkage  is  constrained  by  the  filler  particles,  the 
volume  nominally  occupied  by  the  binder  has  a 
large  void  content.  The  resulting  macropores  in 
the  binder  allow  for  easy  passage  of  evolved  gases 
through  the  binder  matrix  and  account  for  the 
relatively  high  porosity  of  bulk  graphite.  Each 
filler  particle  is  held  in  the  resulting  matrix  by  a 
network  of  binder  bridges  occasionally  interrupted 
by  interconnected  macropores.  The  macropore 
diameter  varies,  but  is  about  the  same  as  the  diam¬ 
eter  of  the  filler  particles.  Aline  of  interconnected 
macropores  may  propagate  a  mechanical  failure 
since  this  is  the  path  of  minimum  cross  section. 

The  microstructure  evidenced  in  Figure  1  exposes 
the  heterogeneous  qualities  of  graphite  in  detail. 
Filler  particles  consist  of  closely  packed  crystal¬ 
lites  and  intercrystalline  micropores.  Significantly, 
the  crystallites  are  not  completely  ordered;  the 
binder  crystallites  are  even  less  ordered,  which 
accounts  for  the  minimal  anisotropy  exhibited  in 
the  thermophysical  properties  of  bulk  graphite. 


Fi|ura1.  Schwiwtie  of  Graphite  Microitructuri 


In  this  paper  ablation  test  data  are  presented 
ior  a  number  of  commercial  graphites  encompass¬ 
ing  fine-  to  coarse-grain  filler  particles.  The 
mean  particle  diameters  varied  between  5  to  100 
microns.  A  typical  particle  diameter  distribution 
determined  by  a  Coulter  counter  is  shown  in 
Figure  2.  For  this  particular  graphite,  50%  of  the 
particles  evidenced  a  particle  diameter  greater 
than  4  microns.  For  the  purpose  of  comparison, 
the  diameter  of  a  grain  of  sand  is  about  100  times 
larger.  The  heterogeneous  character  of  graphite 
is  again  exhibited  in  the  particle  size  distribution. 
Further,  studies  have  been  conducted,  using  X-ray 
diffraction  analysis,  on  the  effect  of  graphitization 
temperature  on  interlaminar  spacing.  Figure  3 
shows  distinctively  different  interlaminar  spacings 
for  the  binder  and  filler.  For  this  particular 
graphite,  the  filler  particles  appear  to  be  com¬ 
pletely  graphitized  at  about  3000 °C;  the  binder 
remains  partially  graphitized.  Complete  graphiti¬ 
zation  occurs  when  the  lattice  spacing  is  about 
3.354A.  Properties  of  the  commercial  graphites 
tested  and  presented  here  are  given  in  Table  I. 
Although  the  density  range  of  the  bulk  graphites 
tested  was  not  large,  appreciable  variations  in 
constituent  material  content,  density,  and  porosity 
are  included.  Nominal  values  cited  are  represent¬ 
ative  of  the  finer  grain  graphites. 


Figutt  2.  Typical  Paiticit  Sizt  Diitributlon 


FiguraS.  Interlaminar  Spacing  ai  a  Function  of 
GrteMtization  Tamparatura  for  Typical  Fina-Grain 
Graphite 


Table  I 

PROPERTIES  OF  BULK  GRAPHITES 


Property  Renge  Nominal  Value 


p  -  Bulk  denelty 
PjT  -  Filler  denelty 
Pq  •  Binder  bulk  denelty 
Kg  -  Maas  fraction  binder 
Cg  -  Binder  poroelty 


1.  65  to 

1.95 

1.8  gm/cc 

Z.  0 

to 

z.z 

2,  IS  gm/cc 

0.8 

to 

1.4 

1.2  gm/cc 

O.Z 

to 

0.4 

0.25 

30 

to 

50 

40% 

Analysis  Model 

Development  of  an  analytical  model  to  repre¬ 
sent  the  complex  and  varied  structure  of  graphite, 
as  schematically  portrayed  in  Figure  1,  is  not 
easily  accomplished.  Rather  convincing  evidence 
has  been  shown  to  destroy  the  concept  of  graphite 
as  a  homogeneous  solid.  It  remains  to  construct  a 
physical  model  whereby  quantitative  information 
concerning  the  surface  erosion  mechanism  of  the 
amorphous  solid  can  be  obtained. To  aid  in  formu- 


’i’Presently,  this  analysis  is  limited  to  loss  of 
particulate  solid  at  the  heated  surface  caused  by 
preferential  mass  loss  of  the  binder  constituent: 
macrospallation  (the  loss  of  larger  scale 
material)  is  omitted,  as  is  very  complex  thermal 
stress  failure. 
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lation  of  a  performance  model,  niany  photomicro* 
(;raphs  of  Kraphitea  exposed  to  hi^h-preaaurc  jet 
heating  conditions  have  been  examined.  Photo¬ 
micrographs  typical  of  the  post-test  condition  of  the 
sectioned  test  models  (subsequently  described)  are 
given  in  Figure  4.  The  cross  section  for  both  a 
coarse-  and  fine-grain  graphite  is  shown  extending 
from  the  heated  surface  to  the  interior.  Indicated 
in  the  figure  is  the  irregular  (roughened)  surface 
consistent  with  the  graphite  mir restructure  shown 
in  Figure  1.  Pockets  at  the  surface  are  evident 
(otherwise  not  seen  in  the  pre-test  models). 
Presumably,  these  are  depressions  caused  by  loss 
of  particulate  solid  at  the  surface.  Particularly 
noticeable  is  the  open  surface  structure  and 
increase  in  macropore  area  in  close  proximity  to 
the  surface.  Attention  should  also  be  given  to  the 
rather  random  distribution  of  the  pores  and  appar¬ 
ent  interconnected  pore  structure. 


0.01  IN.  |.«_ 
COARK  GRAIN 


Figun  4.  Griphitt  Photomkrographi 

The  mechanism  by  which  graphite  surface 
erosion  is  effected  can  then  be  postulated.  At  ele¬ 
vated  surface  temperatures  commensurate  with  the 
high  pressures  addressed  here,  molecular  oxygen 
reacts  chemically  with  the  surface,  independent  of 
the  orientation  of  the  crystallites  (hence,  exposure 
of  edge  planes).  The  extent  of  oxidation  is  deter¬ 
mined  solely  by  the  transport  rate  of  oxygen  to  the 
surface,  and  is  thus  diffusion  controlled.  The 
order  of  the  crystallites,  however,  is  reflected  to 
an  extent  on  the  solid  properties;  entropy  directly 
reflects  the  disorder.  At  surface  temperatures  in 
excess  of  5000°R  (the  threshold  temperature  for 
sublimation)  considerable  quantities  of  carbon  vapor 
are  evolved  from  the  solid.  Sublimation  is  not  con¬ 
strained  to  the  surface:  the  open.  3 -dimensionally 


connected  binder  structure  beneath  the  surface 
affords  a  means  by  which  Internal  sublimation  can 
occur.  Preferential  erosion  of  the  binder  constit¬ 
uent  is  primarily  by  these  two  mechanisms,  expos¬ 
ing  filler  particles  at  the  heated  surface. 


The  phenomenon  causing  loss  of  particulate 
solid  at  the  surface  is  shown  schematically  in 
Figure  5.  Filler  particles  of  mean  characteristic 
dimension,  F,  are  exposed  above  the  receding  sur¬ 
face  at  an  Initial  height,  h.  During  this  time  incre¬ 
ment,  the  particle  recedes  by  an  amount,  Sp, 
caused  from  chemical  consumption  via  chemical 
reactions  with  boundary  layer  air  reactants  and  by 
direct  sublimation.  A  transformed  coordinate  sys¬ 
tem  (y  fixed  in  space)  is  utilized  to  describe 
particle  exposure  in  terms  of  particle  and  binder 
recession  rates.  The  individual  recession  rates  for 
the  binder  and  filler  phases  are,  respectively. 


(1) 

• — 

•E 

M 

(2) 

Figun  S.  Fsrticls  Erosion 


Employing  the  transformed  coordinate  system,  the 
actual  height  of  the  exposed  filler  particle  above 
the  moving  coordinate,  Xg,  is 


Equation  3  provides  definite  information  about 
particle  exposure.  At  this  time,  the  ratio  of  mass 
fluxes  in  Equation  3  is  taken  as  unity,  assuming 
identical  properties  of  the  two  solid  phases  as 
affects  mass  transfer.  Referring  to  Table  I  the 
density  ratio  is  about  0.  5.  Thus,  for  a  particle  not 
acted  upon  dynamically,  and  fully  exposed  (h/F  =  1), 
the  particle  has  receded  half-way  through  the 
exposed  height  and  is  freed  from  the  binder.  In 
the  dynamic  situation,  exposure  of  particles  occurs 
continuously,  and  subsurface  sublimation  causes 
loss  of  binder,  burning  through  the  binder  bridges 
holding  the  particles  as  the  macropores  enlarge. 
During  this  time  the  particles  exposed  have  partially 
receded.  Finally,  the  filler  particles  are  restrained 
only  by  geometric  interlocking  with  adjacent 
particles,  and  presumably,  would  be  blown  free  by 
moderate  pressure  or  shear  forces.  The  erosion 
of  individual  filler  particles  is  consistent  with 
available  ablation  data,  since  the  particles  are  too 
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■mall  (0.  02  to  5  microna)  to  be  obaerved  by  normal 
photographic  techniquea,  and  the  eurface  would 
appear  to  ablate  uniformly.  Preferential  loae  of 
binder  and  expuleion  of  filler  particle*  from  the 
iurface  has  been  obeerved  by  Wachi  in  free 
vaporization  atudiea  conducted  on  a  number  of 
graphite  teat  aamplea. 

The  remainder  of  thia  paper  ia  devoted  to 
determination  of  the  fraction  of  particulate  aolid 
eroded  from  the  receding  aurface  in  terma  of 
environmental  exposure  conditions,  Tn  accomplish 
identification  of  the  thermomechanical  recession 
component  from  surface  erosion  measurements  of 
the  test  data  we  first  describe  the  thermochemical 
recession  in  some  detail.  The  probable  uncer¬ 
tainty  in  thermodynamic  data  employed  on  surface 
recession  rates  is  evaluated  in  a  description  of 
sublimation  rates.  Loss  of  particulate  solid  at  the 
heated  surface  is  then  addressed  as  affects  mass 
and  energy  transfer.  Correlation  of  the  experi¬ 
mental  data  in  terms  of  thermochemical  erosion 
enables  quantitative  identification  of  the  solid 
fraction  eroded. 

III.  Carbon  Thermodynamics 


Considerable  interest  and  research  have  been 
focused  on  the  identification  and  thermodynamic 
characterization  of  high-temperature  carbon  vapor 
species,  and  elucidation  of  the  carbon  vapor  equa¬ 
tion  of  state.  Up  to  now,  the  thermodynamic 
properties  of  all  polyatomic  carbon  vapor  species 
larger  than  Cj  have  been  calculated( ^2,  13,  l4)  from 
the  (L-F)  "  "  increments  of  Pitzer  and  Clementi(12), 
which  were  fitted  to  ^eld  a  triple-point  pressure  of 
100  atm  at  4000‘’Kn5)  LU^ewise,  the  equation-of- 
state  calculations'”’  have  failed  to 

include  the  heavier  molecular  weight  species 
Ci7  to  C30,  despite  evidence  of  their  existence^' 
and  likelihood  of  their  prominence  at  high  pressures 
and  temperatures.  The  carbon  ring  molecules  have 
also  been  excluded.  In  this  paper,  both  linear  and 
polygonal  ring  structures  of  are  considered  up 
to  C30.  The  resonance  energies  and,  ultimately, 
the  thermodynamic  properties  of  the  molecules  are 
calculated  by  theoretical  methods  substantially 
different  from  those  used  by  Pitzer  and  dementi'*^' 
and  later  by  Strickler  and  Pitzer^^^^  The  theoret¬ 
ical  procedures  are  internally  self-consistent,  and 
yield  results  which  are  in  good  agreement  with  the 
observed  sublimation  data  and  recent  triple -point 
temperature  measurements  of  graphite. 

Vapor  Properties 

The  molecules  studies  were  linear  chains  with 
4  to  30  carbon  atoms,  and  regular  polygonal  rings 
with  3  to  30  carbon  atoms.  Figure  6  depicts  the 
assumed  valence  bond  configuration  for  each  type 
of  structure.  Bonding  is  assumed  to  consist  of 
localized  <s  bonds  and  delocalized  tt  bonds  construc¬ 
ted  from  linear  combinations  of  2p  atomic 
orbitals.  The  carbon-carbon  distance  is  taken  as 
1.  28A:  the  bond  length  in  03'^^'  and  the  C=C  dis¬ 
tance  in  carbon  suboxide,  C3O2  (2^).  The  tt 
orbitals  in  the  linear  chains  are  delocalized  with 
respect  to  the  xz  and  yz  planes,  with  the  2p^ 
atomic  orbitals  pointing  in  the  x  and  y  directions. 
The  ring  compounds  have  a  set  of  ir  orbitals 
formed  from  2p  atomic  orbitals  oriented 
perpendicular  tS  the  plane  of  the  ring  and 
another  set  formed  from  2p^  orbitals  oriented 
radially  in  the  plane  of  the  ring.  Except  in  large 


rings,  the  latter  do  not  exhibit  true  n-bund 
behavior  because  of  increasing  a  and  p  orbital 
interaction  with  increased  bending  of  the  carbon- 
carbon  chain,  I  The  electronic  ground^state  of 

linear  C  ia  /  if  n  ia  odd,  and  /  „  if  n  is 
n  ^g  K 

even.  In  the  ring  molecules,  the  singlet  state 

occurs  where  n  is  even,  and  the  triplet  state 

where  n  is  odd. 


FigunB.  Vilancs  Bond  Structure  for  Linear  Chain  and 
C^oPolygonal  Ring 

The  heats  of  formation  of  the  polyatomic 
carbon  species  were  calculated  from  the  relations 
of  Strickler  and  Pitzer(21). 

Linear  Molecules 

"  ^“c3(L)  ^  ^C^(L)  ■  ^CjlR) 

+  (n-3)E.  -  23/n  {1  -  Mod(n,  2)}  (4) 


Ring  Molecules 


AH®  =  nE.  +  R^  -  23/n  Mod(n.2)  (5) 

n  n' 

The  operator  Mod  (n, 2)  in  Equations  1  and  2  is 
defined  as  n  -  2  j  n/2}  and  is  a  correction  for  the 
difference  in  energy  between  the  triplet  electronic 
ground  state  and  the  slightly  higher  singlet  state 
associated  with  the  calculation, 

Ej  (the  energy  per  interior  carbon  atom)  was 
set  equal  to  43.  4  kcal/mole,  the  value  calculated 
by  Pitzer  and  Clementi(l2)  from  the  energy  differ¬ 
ence  in  allene-ethylene,  corrected  by  the  (L-^) 
method  for  the  effect  of  hyperconjugation. 
was  computed  to  be  202.  6  kcal/mole  from 
Drowart' 8(24)  rnass  spectrometric  data,  the  Third 
L,aw  relationship,  and  the  thermodynamic  proper¬ 
ties  of  linear  C3  calculated  here,  A  good  method 
is  not  available  for  estimating  the  energy  of  the 
terminal  at^ms  in  a  linear  chain;  the  observed 
value  of  AHq^(L)  is  taken  as  the  reference  point. 

The  resonance  energy,  R,  of  an  unsaturated 
molecule  is  defined  as  the  difference  in  electronic 
energy  of  the  actual  delocalized  structure  and  a 
hypothetical  localized  double-bond  configuration. 
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The  eneriiiei  of  the  delocalised  and  localised  mole¬ 
cular  urbitala  were  calculated  by  the  extended 
llUckel  technique  of  Hoffmannt^S),  in  particulart 
the  delocalised  energy  waa  obtained  by  solving  the 
cun'plete  set  of  HUckel  equations 


n 


i  1 


j  1.2,...n  (6) 


using  as  a  basis  set  n  2s  and  3n  2p  carbon  Slater 
orbitals  with  exponents  of  1,  625.  The  matrix  ele¬ 
ments,  Hij,  were  evaluated  by  the  relations  in 
Reference  25,  the  overlap  integrals,  Sjj,  were  cal¬ 
culated  directly.  The  energy  for  the  hypothetical 
localised  structure  was  obtained  from  Equation  6 
by  setting  some  of  the  Hij  and  Sij  matrix  elements 
of  the  2p^  atomic  orbitals  equal  to  zero,  so  as  to 
have  a  localized  structure, 


Because  of  the  extremely  large  computer  run 
times  required  for  the  larger  species  (e.  g.  ,  1.3 
hours  for  C24),  the  resonance  energies  of  the 
large  molecules  were  obtained  by  the  interpolation 
or  extrapolation  of  smooth  curves  constructed  from 
the  calculated  resonance  energy  data  of  linear  and 
ring  C3  to  Cjg  and  C24.  Small  errors  introduced 
by  this  procedure  are  of  no  consequence  except, 
perhaps,  at  very  high  temperatures  and  pressures 
outside  the  range  of  this  study.  The  delocalized 
energies  of  several  of  the  carbon  species  were 
taken  directly  from  Hoffmann's  orbital  energy 
data(2°).  These  energies  were  obtained  for  a 
C  C  distance  of  1.  30A  for  the  C_  species;  however, 
the  use  of  I.  30A  instead  of  1. 28A  has  no  signifi¬ 
cant  effect  on  the  calculated  energy.  Table  II  lists 
the  resulting  heats  of  formation  of  both  the  linear 
and  ring  carbon  vapor  species;  the  values  of 
for  Cj  and  €3  were  obtained  from  Reference  24. 
The  calculated  heats  of  formation  for  cyclic  C3  and 
C4  were  found  to  be  too  low  due,  presumably,  to 
an  underevaluation  of  the  strain  energy.  These 
molecules  are  not  included  in  these  calculations. 

In  addition,  another  set  of  AHp  values  was  calcu¬ 
lated  using  the  method  of  Striclder  and  Pitzer^^^K 
(A  revised  version  of  the  original  Pitzer  and 
Clementi  calculation,  )  These  alternate  calcula¬ 
tions  were  performed  with  a  strain  energy  for  the 
ring  compounds  of  154.  5/n  kcal/mole  based  on  a 
fundamental  bending  frequency  of  63  cm"^  for 
C3O2,  and  a  heat  of  formation  of  202.  6  kcal/mole 
for  linear  C3.  The  AHp's  obtained  by  the  Hoffman 
method^^^)  were  as  much  as  1  507#  larger  for  cyclic 
Cg,  and  an  average  of  16%  and  18%  larger  for 
cyclic  and  linear  C20  lo  Cjn,  than  the  values 
obtained  by  Strickler  and  Pitzer.  These  differ¬ 
ences  may  be  partially  accounted  for  by  the 
different  definitions  of  H^:  utilized  in  the  two 
methods.  However,  the  extended  IRickcl  theory  of 
Hoffman,  unlike  the  LCAO-MO  method  of  Strickler 
and  Pitzer,  treats  the  entire  secular  determinant 
(including  non-nearest  neighbor  interactions),  and 
takes  into  account  the  precise  atomic  coordinates 
of  the  various  atoms.  Other  considerations  being 
equal,  the  above  should  result  in  orbital  energies 
which  are  more  accurate.  This  is  particularly 
true  in  the  case  of  the  intermediate-size  ring 
s  pe  c  i  e  s . 

Thermodynamic  properties  of  the  carbon  vapor 
species  were  calculated  by  the  well-known  relations 
of  statistical  mechanics.  The  electronic  states  of 
Cj  were  obtained  from  Reference  27.  The  spectro¬ 
scopic  constants  of  €3  were  taken  from  References 


28  and  29.  The  conetante  for  C3  war*  obtained 
from  Reference  30;  a  bending  force  conetant  of 
63  cm'l  waa  employed  in  the  calculation.  The 
spectroscopic  constants  (i.e.,  moments  of  inertia 
and  fundamental  vibrational  frequenciea)  of  all 
other  carbon  vapor  species  were  calculated  and 
derived  directly.  Equations  for  calculating  the 
fundamental  vibrational  frequencies  of  the  linear 
and  ring  molecules  were  derived  from  a  normal 
coordinate  analysis  of  the  molecule.  Assumptions 
were  made  in  the  derivation  that  the  valence  force - 
field  consisted  of  bending,  stretching,  and  (in  the 
case  of  rings)  torsional  motion;  and  that  all 
motions  wore  harmonic  in  nature,  and  mutually 
independent.  The  equations  for  the  fundamental 
frequencies,  (cm*l),  associated  with  each  type  of 
structure  and  molecular  motion  are  given  below; 

Linear  Molecule 

(a)  Stretching  motion 

{2itcu)^  =  (4Xg/m)  sin^  (k/2)  (7) 

(b)  Bending  motion 

(2itcw)^  =  (l6Xg/m)  sin^  (k/2)  (8) 


1/2  1,1 


^11  '  (8Xg/m)(r/s)^  sin"*  (n/nXl  +  cos  k)^ 

+  (2Xg/m)(r/B)^  sin^  (2iT/n)(3  +  cos  2k  -  4  cos  k) 

^22  ’  (8^/”')(r/B)^  Bin^  (n/r)  cob^  (TT/n)(  1  -  cob  k)  . 

+  (2Xg/m)(r/B)^  (2  sin^  n/n)^  ( 1  -  cos  2  k) 

f  =  (8Xg/m)(r/8)^  Bin^  (  n/n)  cob  ( Tr/n)  sink 

+  (2X_/m)(r/B)^  8in(2TT/n)(2  sin^ir/n)^ 

15 

2  2  2 

(Bin2k-8ink)  (2ncw)  =  (8X.j,/mB  sin  2TT/n) 

2 

x(  1  -  COB  k)(coB  k  -  cos  2  ir/n)  (10) 


^S’  and  Xj  are  the  respective  stretching,  bend¬ 

ing  and  torsional  force  constants  for  both  chain  and 
ring  structures.  The  bending  mode  frequencies  in 
Equation  8  are  doubly  degenerate,  and  Equations 
9  and  10  are  zero  for  three  distinct  values  of  j. 
Consequently,  Equations  7  through  10  yield  a  total 
of  2n-5  and  2n-6  values  (not  necessarily  distinct) 
of  u)  in  the  case  of  the  linear  and  ring  molecules, 
respectively,  in  accordance  with  the  principle  of 
equi partition.  The  force  constants  Xg  and  X3  were 
assigned  the  values  11.83  x  10^  dyne  cm'l  and 
5,  100  dyne  cm*l.  These  are  the  magnitudes  of  the 
carbon-carbon  stretching  force  con8tant(32)  and 
the  central  carbon  bending  force  constant^^^)  of 
C3O2  •  X J  was  set  equal  to  3.36x  10"12  erg/rad, 
the  value  of  the  torsional  force  constant  of 
benzene'^^K 
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Tabu  II 

PROPERTIES  OF  CARBON  SPECIES 
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294.0 
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395,  1 

245.9 

178.0 

113.3 
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49.4 

486.4 

326.3 

254.9 
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442.! 

357.4 

210.0 

44.9 

911.4 

338.8 

242.  1 

188.9 

Cu 

474.4 

370.5 

212.2 

54.2 
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327.  5 

252.5 
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‘=19 
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230.2 
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595.0 
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629.) 

431.  5 

347.  9 

268.8 
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445. 1 

257.1 

71.8 
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944.7 
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94.2 
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‘=30 
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Solid-Vapor  Equilibrium 

Table  II  lists  the  standard -state  free  energies 
of  formation  of  the  individual  carbon  species  at 
three  different  temperatures,  calculated  from 

(UI 


and  should  be  included  in  calculations  performed 
above  4000*K.  Except  for  C3,  and  possible  C5, 
the  odd-even  alternation  in  stability  computed  by 
Pitser  and  dementi,  and  observed  by  Baun  and 
Fischer  up  to  C^,  is  not  found  in  the  linear  chain 
data  of  Figure  7.  However,  the  measurements  of 
Baun  and  Fischer  were  made  on  ions,  and  the  trend 
in  stability  may  shift  when  applied  to  the  neutral 
counterparts.  These  investigators  obtained 
increased  concentrations  of  doubly  charged  ions 
C^'*' with  n  =  11,  15,  19  (whose  structures  are 
believed  to  resemble  the  closed-shell  configura¬ 
tions  of  the  neutral  cyclic  species  Cjo,  C14,  C]g) 
which  have  superior  stabilities  (as  shown  in 
Figure  8).  In  Figure  7,  an  even-odd  alternation  is 
evident  above  Cio  at  4000 ‘K;  the  effect  appears  to 
propagate  to  smaller  values  of  n  with  increasing 
temperature. 


Figure  8.  Pirtiil  Pretwrei of  Ring  Cirbon  Vipor  SpKitt 


related  to  the  partial  pressure,  P  ,  of  C  by  the 
equilibrium  constant  ^n 


RT 


InKp 


(12) 


and 


In  Figure  9,  the  total  pressure  of  the  vapor  in 
equilibrium  with  graphite 


P  ^ 


n=  1 


(L)  ^  *C  (R) 

n  n 


(14) 


Kp  :  (13) 

The  free  energy  of  graphite  solid  used  in  Equation  11 
was  taken  from  JANAF^^®). 

Figures  7  and  8  depict  the  logarithms  of  the 
partial  pressures  of  the  respective  linear  and  ring 
molecules  in  equilibrium  with  graphite.  The  plots 
show  linear  C3  is  the  most  stable  species  and 
linear  C5  somewhat  more  stable  than  and  Cg,  in 
agreement  with  previous  mass  spectrometric  find- 
ings(*°*  It  is  also  apparent  that  the  contribu¬ 

tion  of  ring  molecules  is  negligible,  and  that  larger 
linear  molecules  become  very  important  at  4500 *K, 


is  shown  and  compared  to  that  calculated  by  Pitzer 
and  dementi,  and  Duff  and  Bauer.  It  is  apparent 
the  data  of  Pitzer  and  Duff  do  not  agree  well  with 
the  experimentally  determined  sublimation  point  of 
Thorn  and  Winslow^^®'  and  the  recent  triple -point 
temperatures  of  Schoessow^^")  and  Fateeva'^'). 

The  results  obtained  match  the  sublimation  point 
and  lie  between  the  more  recent  triple -point  tem¬ 
peratures.  This  calculation  method  is  advanta¬ 
geous  because  (1)  it  generates  the  thermodynamic 
properties  of  both  linear  and  ring  structures  in  a 
self-consistent  manner,  and  (2)  the  resulting  data, 
used  in  conjunction  with  the  new  heat  of  formation 
values,  are  consistent  with  the  latest  experimental 
data  without  recourse  to  alteration. 
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Figura  9.  Graphiti  Equilibrium  Vapor  Pmsura 


Because  the  free  energies  of  formation  ulti¬ 
mately  control  the  state  properties  of  the  system, 
the  influence  of  AFp  on  the  state  properties  was 
stiuliecl.  Calculation  of  these  properties  was 
repeated  with  the  free  energy  of  each  species  with 
n  >  3  perturbed  by  ±20%.  The  limits  were  adopted 
on  consideration  of  the  probable  cumulative  uncer¬ 
tainty  in  the  calculated  heats  of  formation  and  free 
energyv.  The  uncertainty  limits  reflect  a  study  of 
the  influence  of  key  spectral  quantities  andmethods 
on  these  properties.  As  discussed  previously,  our 
heats  of  formation  for  the  large  molecules  are 
about  18%  greater  than  the  Strickle r  and  Pitzer 
values,  accordingly  increasing  AFp.  In  calcula¬ 
tion  of  the  free  energy,  an  evaluation  was  made  of 
the  effect  of  the  bending  force  constant  (^p)  in 
perturbing  the  results.  The  bending  force  con¬ 
stant  contributes  more  strongly  to  the  theimody- 
namic  properties  than  the  stretching  and  torsional 
force  constants.  Figures  10  and  11  show  the  influ¬ 
ence  of  on  the  free  energy  foe  selected  linear 
and  cyclic  molecules.  The  free  energy  and  hence 
the  free  energy  of  formation  decrease  markedly  as 
the  bending  force  constant  is  decreased  below 
‘i,  100  dyne  cm"^.  A  50%  decrease  in  Xg,  for 
example,  lowers  the  free  energy  of  formation 
approximately  lO'!',,  at  3000  °K  and  20%  at  4000  "K 
for  the  larger  species.  Thermodynamic  calcula¬ 
tions  made  here  used  the  bending  force  constant 
for  C3O2  (5,  100  dyne  cm"l).  This  molecule  was 
selected  because  of  the  similarity  in  bonding  with 
and  because  the  carbon  atoms  are  less  apt  to 
exhibit  terminal  end  effects  than  the  atoms  in  Cj. 

Figures  9,  12  and  13  show  the  effect  of  the 
uncertainty  in  AF^  on  the  respective  equilibrium 
vapor  pressure,  average  molecular  weight  of  the 
vapor,  and  sublimation  energy  of  graphite.  The 


The  much  smaller  uncertainty  in  the  properties 
for  Ci|-Cj  is  omitted  in  these  calculations. 


total  equilibrium  pressure  (Figure  9)  (hence  the 
number  of  moles  of  each  species)  increases  sub¬ 
stantially  as  the  free  energy  of  formation  is 
decreased.  The  effect  on  the  average  molecular 
weight  and  sublimation  energy  is  likewise  pro¬ 
nounced.  _The  trend  in  Figure  12  is  for  larger 
values  of  Mwith  decreasing  AF|<.  This  is  due  to 
the  larger  increase  in  population  of  the  heavier 
molecular  weight  species  in  relation  to  the  lighter 
species.  In  Figure  13,  the  sublimation  energy  of 
graphite,  which  is  defined  as 


AH 
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^  (L) 

Z—  »C(L) 
n=  1  n 


^Cn(R) 

P  (R)  ■ 

n  s 

tends  to  decrease  significantly  as  the  free  energy 
of  formation  is  increased. 


Figura  10.  Influinci  of  Bonding  Forci  Conitint  on  Frai 
Enirgy  it  3,000*’K 
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IV.  Mage  and  Energy  Transfer 

Application  of  the  governing  conservation 
equations  to  an  ablating  body  with  internal  gas 
generation  has  been  well  studied.  Basically,  the 
differential  equations  used  to  describe  thermo¬ 
chemical  erosion  in  multicomponent-chemical 
reacting  flows,  >tem  from  the  hypersonic  approxi¬ 
mations  of  Lees(^°^  formulations  of  Scala(j/,  and 
those  derived  by  Dorrance(39)^  When  similarity 
constraints  are  retained  (e.g.,  at  the  stagnation 
point,  and  for  certain  body  geometries)  the  differ¬ 
ential  equations  reduce  to  differences  only  in  the 
thermodynamic  and  transport  properties  employed 
and  the  numerical  techniques  used.  Lees  adopts 
the  effective  binary  diffusion  coefficient  model  with 
unity  Lewis -Semenov  number,  obtaining  a  direct 
solution  for  mass  and  energy  transfer.  Seals 
obtains  numerical  solutions  to  the  differential 
equations  employing  transport  and  thermodynamic 
properties  for  each  chemical  species  included  in 
his  reaction  system,  which  appears  limited  by 
present-day  computer  storage  capacity.  An 
approximation  to  large  multicomponent  reaction 
systems  has  been  developed  by  Kendall,  et  alH^^ 
utilizing  the  differential  equations  of  Dorrance, 
and  achieving  more  versatility  than  either  Lees  or 
Seals.  It  remains  that  significant  differences 
between  the  several  sets  of  solutions,  when  using 
identical  thermodynamic  properties,  has  yet  to  be 
demonstrated.  Thus,  the  numerical  results 
calculated  here  retain  the  similarity  approxi¬ 
mations  of  Lees. 

Mass  Transfer 


The  elemental  mass  balance  of  a  material 
consisting  of  i-chemical  elements  whose  surface 
recedes  at  a  rate  s  and  undergoing  loss  of  solid 
material,  nix',  at  the  surface  (Figure  14)  is  given 
by 

m  K.  +  m  R.  -  ml'  R.*  -  (pv).  =  0  (16) 

g  1  SI.  i  1 

g  s  w 

where  mg  Kj.  and  mg  Kig  are  the  mass  fluxes  of 
the  element  rinto  a  control  volume  fixed  at  the 
receding  surface  owing  to  an  internal  gas  source 
and  solid  material  transfer,  respectively.  Mass 
leaves  the  control  volume  normal  to  the  wall  by  a 
net  flux  of  gas,  (Pv)i  ,  and  by  loss  of  solid, 
rh*  Ri*.  The  mass  flux  of  the  i*^  element  normal 
to  the  wall  accrues  by  a  summation  of  individual 
flux  components  consisting  of  the  internal  gas 
source,  surface  chemical  reactions  with  the 


■(‘Comparisons  between  the  equal  and  unequal 
diffusion  coefficient  model  of  Lees  and  Scala  has 
been  made  by  Hearne,  etaL.WD  and  Dolton, 

etal.(14). 


external  reactive  gas  (boundary  layer  and  internal 
gas  species)  and  by  direct  sublimation.  For 
graphite  in  an  air  boundary  layer  the  following 
conservation  statements  are  rnade  at  the  wall, 

(c.f. ,  I..ees,  Kratsch,  etal.'^O. 


Rj^  +  Kq  +  =  1  (17) 

WWW 


(pv)  .  pv  =  pv^R 

^W  w 

(pv)  =0  =  pv  R,.. 

w  0^ 

(pv)j^  -  0  -  pv^R^ 

w  w 


The  gradient  of  the  individual  chemical 
elements  at  the  wall  is  given  by,  e.  g. ,  carbon 


and  for  unity  Lewis-Semenov  number  the  convec¬ 
tive  heat  transfer  to  the  wall  by  conduction  and 
diffusion  is,  upon  introduction  of  Stanton  number 


Thus,  Equation  16  can  be  solved  by  introducing 
Equations  18,  21,  and  22  and,  after  rearrangement, 
the  elemental  mass  fraction  of  carbon  at  the  wall 
is  given  by 


rh^  K„  m  K„ 

B  ^  ®  ^8 


K, 


p  u  C„ 
e  e  e  H 


P  u  C„ 
e  e  H 


(23) 


1  + 


(pv) 

w 

P  u  C„ 
e  e  H 
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where  for  graphite  in  an  air  boundary  layer 


ni  +  ni 
g  e 


omit  carbon  molecules  larger  than  C3.  Upon  defi¬ 
nition  of  a  surface-energy  balance  and  the  in-depth 
heat  conduction,  sublimation  rates  are  directly 
coupled  to  surface  temperature  and  pressure  and 
the  local  aerodynarhic  heating. 


0.  768 


--  0.232 


=  1 


The  fraction  of  solid  entering  the  control 
volume  which  is  lost  at  the  surface  as  particulate 
solid  is  defined  from  Equation  24 


m-'  {  m 

8  S 

Inspection  of  Equation  23  shows  that  the  mass 
fraction  of  the  element  carbon  at  the  wall  is  deter¬ 
mined  solely  by  the  mass  flux  of  elemental  carbon 
injected  at  the  surface  as  vapor,  reflected  by  the 
difference  in  fluxes  entering  and  leaving  the  con¬ 
trol  volume.  Equations  similar  to  Equation  23  for 
the  element  oxygen  and  nitrogen  provide  the  neces¬ 
sary  relations  satisfying  Equation  17  in  terms  of 
the  individual  fluxes  depicted  in  Figure  14.  A 
stoichiometric  statement  then  directly  relates  the 
elemental  mass  fractions  to  the  chemical  composi¬ 
tion  of  the  combustion- sublimation  species  at  the 
ablating  surface 

^  V.  M. 

•^i  I  -4r  ^ 

j  •> 

Numerical  solutions  to  the  governing  equations 
were  made  invoking  the  chemical  equilibrium  con¬ 
straint.  The  chemical  reaction  system  included 
sublimation  species  for  the  ring  and  chain  mole¬ 
cules  given  in  Table  II,  and  in  addition  air  species, 
and  carbon-air  reaction  products  listed  below* 


SURFACE  lEMPERATURE  TRI 

Figuri  16.  Gnphitt  Sublimition  Riti 
Transfer 


Energy  transfer  to  a  surface  undergoing  mate¬ 
rial  consumption  via  thermochemical  erosion  and 
thermomechanical  loss  of  solid  is  given  by  an 
energy  balance.  Figure  I6,  The  heat  conducted 
by  the  solid  away  from  the  surface  reflects  the 
energy  fluxes  leaving  and  entering  the  control 
volume  expressed  in  the  following: 

*^conv  ^rad  *  ^rad  *  ^w  ~  " 

in  out 


where 


-  rh 

^cond  8  K  8  8 


2«j«i 

j 

H.=  JCpdT+^H° 


^2^2 

C4N2 


^conv  ^rad  ^rad 

in  out 


Graphite  sublimation  rates  over  a  large  pres¬ 
sure  range  are  shown  in  Figure  15  normalized  to 
the  familar  diffusion- controlled  oxidation  rate 
(riiD  =  0.  174  PgUg  C|-{).  In  the  first  several 
decades  of  pressure  the  threshold  sublimation 
temperature  increases  with  increasing  pressure 
reflecting  the  carbon- vapor-pressure  diagram 
shown  in  Figure  9.  The  10  ^  to  10^  atm  pressure 
range  is  of  special  interest;  the  threshold  sublima¬ 
tion  temperature  is  compressed  due  to  sublimation 
of  large  cjuantities  of  polyatomic  carbon-vapor 
species.  Thus,  the  high-pressure  sublimation 
regime  of  graphite  occurs  over  a  fairly  narrow 
surface  temperature  range  (circa  7, 000- 8,000° R)in 
contrast  to  Scala's(3)  results  (7, 000- 1 0, 000  °  R)which 


‘>cond  l^gHg  hsHj 

Figure  16.  Energy  Trentfer 


The  summation  in  Equation  28  is  made  over 
all  species  in  the  gas  at  the  wall  with  concentra¬ 
tions  defined  by  the  mass-transfer  rates  and  appli¬ 
cation  of  the  equilibrium  constant.  For  the  case  of 
Le  =  1  the  net  heat  transfer  to  the  surface  may  be 
expressed  as  (at  this  time  rhg  is  set  equal  to  zero) 

q  J  =  p  u  C„  He  +  /l-f  \  B'  H  -  /l+B'  \H 
^cond  '^eeH  Sg  I  si  ss  I  wf  w 


Thermodynamic  property  data  for  these  species 
were  taken  from  the  JANAF  tables. 


+  ^rad  -  '^"^w 
in 


10 


when  the  usual  definition  of  the  mass  transfer 
parameter  is  made 


B'  = 


m 

p  u  C„ 
e  H 


(30) 


The  reduction  in  Stanton  number  due  to  mass 
injection  is  functionally  dependent  on  blowing  rate, 
condition  of  the  boundary  layer,  and  properties  of 
the  injected  species 


C 


r  ='(»'»■  »i„j)  <’» 


Equation  29,  when  coupled  to  a  standard  finite 
difference  solution  to  the  conduction  equation 
enables  a  complete  description  of  the  surface 
erosion  of  graphite  including  solid  loss.  It  is 
observed  that  the  effect  of  solid  removal  on  the 
surface  energy  balance  is  to  reduce  the  net  energy 
transfer.  For  example,  at  a  surface  temperature 
of  7,  500*R  the  B'gHg  term  is  about  1, 700  Btu/lb 
at  100  atm  pressure;  Equation  29  is  reduced  only 
by  850  Btu/lb  for  fg  =  0,  5,  In  comparison,  the 
enthalpy  of  the  combustion- sublimation  products  at 
the  ablating  surface  is  about  4,000  Btu/lb  as  shown 
in  Figure  17.  The  enthalpy  of  the  gaseous  products 
(given  by  Equation  28)  evidences  the  exothermic  CO 
chemical  reaction  product  below  S,000°R  surface 
temperature.  The  increase  in  enthalpy  above 
5,000“R  is  coincident  with  sublimation;  highly 
endothermic  carbon  vapor  species  evolve  at  the 
surface. 


Figurt  17.  Enthilpy  of  Combuition-SuElimition  Products  it 
Abliting  Gnphiti  Surfict 


Figura  18.  Chitnicil  Composition  it  Griphiti  Surfici 

It  is  interesting  to  examine  how  the  ablated 
mass  of  carbon  is  partitioned  among  the  chemical 
composition  shown  in  Figure  18.  In  Figure  19  the 
elemental  mass  fraction  of  carbon  is  shown  extend¬ 
ing  from  the  diffusion- controlled  oxidation  regime 
well  into  the  sublimation  regime.  As  indicated, 
the  constant  mass  fraction  of  0.  15  reflects  the 
carbon  mass  in  CO.  At  about  6,500°R  nitrogen 
reacts  forming  the  cyano  radical  and  cyanogen; 
only  a  small  mass  of  carbon  is  ablated  because  of 
these  reaction  products.  At  about  the  same  tem¬ 
perature  (6,500°R)  contributions  due  to  Cj  to  C3 
occur,  increasing  the  mass  fraction  of  ablated 
carbon.  With  increased  temperature  larger  carbon 
vapor  species  sublime  from  the  surface  appearing 
first  as  species  having  from  4  to  10  carbon  atoms 
and  then  in  species  having  from  11  to  20  carbon 
atoms.  For  the  100-atm  pressure  condition 
shown,  little  contribution  from  vapor  species  with 
greater  than  20  carbon  atoms  occurs.  More 
clearly,  than  indicated  in  Figure  18,  it  is  apparent 
that  the  additive  small  concentrations  of  the  poly¬ 
atomic  carbon  vapor  species  is  quite  appreciable. 


Detailed  results  of  the  high-pressure  sublima¬ 
tion  regime  are  shown  in  Figure  18.  The  chemical 
composition  at  the  ablating  graphite  surface  is 
given  at  100  atm  pressure.  In  addition  to  the  air- 
carbon  reaction  products,  a  complex  set  of  subli¬ 
mation  species  is  evidenced.  The  triatomic 
carbon-molecule  is  most  abundent;  however,  a 
large  number  of  polyatomic  molecules  of  rather 
low  concentration  are  seen.  Other  than  C3  the 
prevalence  of  any  single  carbon  vapor  species  is 
not  evident.  Although  vapor  species  having  ring 
structures  were  included  in  the  numerical  compu¬ 
tations,  the  concentrations  shown  are  due  only  to 
linear  structured  molecules.  Conclusively,  the 
ring  molecules  can  be  safely  neglected  in  perform¬ 
ance  analyses  for  graphitic  materials.  Omission 
of  the  linear  carbon  species  at  high  pressures 
examined  here  may  be  quite  critical  as  affects 
mass  transfer. 


Figura  19.  Pwtition  of  ElinNntal  Mm  Fraction  of  Carbon 
with  Surfaci  Tamporatura 

Sensitivity  Analysis 

The  uncertainty  limits  established  for  the 
carbon-vapor  species  as  affects  mass  transfer  and 
energy  transfer  can  be  assessed  employing  the 
conservation  equations  previously  described.  Fig¬ 
ure  20  shows  the  normalized  sublimation  rate  for 
several  pressures  including  ±20%  uncertainty  in 
free  energy  of  formation  for  the  carbon-vapor 
species.  The  uncertainty  band  reflected  in  Fig¬ 
ure  20  is  consistent  with  the  vapor  pressure  dia¬ 
gram  shown  in  Figure  9;  the  higher  pressure 
evidences  a  greater  uncertainty  band  due  to  the 
greater  concentration  of  sublimation  species.  For 
comparative  purposes,  results  calculated  by 
Scala(^)  are  shown  for  a  normalized  sublimation 
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rate  of  4.0.  Also  shown,  are  results  by  Dolton 
et  al.  at  1  atm  pressure.  He  includes  vapor 

species  up  to  employing  empirically  deter¬ 

mined  free  energies  forced  to  fit  a  presumed 
4,  000°K  tripli'-point  temperature.  It  is  seen  that, 
at  100  atm  pressure  and  for  constant  surface 
temperature,  sublimation  rate  can  be  as  much  as 
3  times  larger  than  the  nominal  value.  The  effect 
of  uncertainty  in  thermodynamic  properties,  dim¬ 
inishes  with  decreased  pressure,  owing  to  the 
smaller  vapor-pressure  contribution  of  the  high 
molecular  weight  vapor  products.  Corresponding 
results  showing  the  effect  of  the  combustion- 
sublimation  products  on  wall  enthalpy,  are  shown 
in  Figure  21.  The  composite  of  these  results 
enable  an  assessment  of  the  uncertainty  on  surface 
recession  rate  to  be  made. 


Figuri  20.  Graphiti  Sublimation  Rata  Saniitivity  to  i20K 
Uncanainty  inApf 
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Figure  21.  Santitivitv  of  Enthalpy  of  CombuHion-Sublimation 
Producti  It  Ablating  Grephitt  Surfica  to  ±20% 
Undartainty  in  AFp 

The  calculated  surface  recession  rate  for 
environmental  conditions  we  have  studied  in  cor¬ 
relation  of  high  pressure  test  data  is  shown  in  Fig¬ 
ure  22. 


Figure  22.  Tharmochamical  SurfKa  Raemion  Rata 


At  the  highest  stagnation  enthaipy  shown 
(10,000Btu/lb)  caiculated  surface  temperature  is 
T.SOO'R  at  100  atm  pressure.  The  sensitivity  of 
surface  recession,  owing  to  the  uncertainty  band 
on  properties,  is  small.  This  sensitivity 
decreases  with  both  reduction  in  stagnation  enthalpy 
and  pressure,  evidencing  the  effect  shown  in  Fig¬ 
ure  20  on  sublimation  rate.  In  the  numericai 
results,  we  have  obtained  both  transient  and  steady- 
state  results,  including  loss  of  solid  at  the  heated 
surface.  In  the  transient  caicuiations  it  was  noted 
that  the  surface  temperature  approached  a  steady- 
state  value  in  less  than  0.  1  second.  For  the  most 
part,  steady- state  predictions  were  made  to  cor¬ 
relate  the  data  because  of  the  large  numerical 
results  required.  Calculated  surface  temperatures 
including  the  effect  of  solid  eroded  is  discussed 
subsequently  in  correlation  of  the  arc  jet  data. 

Even  though  a  large  uncertainty  in  the  thermo¬ 
dynamic  properties  of  the  carbon  vapor  species 
was  included  (resulting  in  almost  a  decade  uncer¬ 
tainty  in  carbon-vapor  pressure)  the  effect  on 
surface  recession  is  less  than  10%  within  the 
environmental  conditions  discussed  in  this  paper. 
This  result  occurs  because  of  the  counter-balancing 
effect  of  increased  endothermic  heat  absorption 
when  greater  quantities  of  carbon  vapor  are 
evolved  at  the  surface.  Thus,  when  using  a  self- 
consistent  model  for  calculating  the  vapor  species 
thermodynamic  properties  with  associated  uncer¬ 
tainty,  the  sensitivity  on  surface  recession  is 
comparable  to  the  accuracy  for  which  local 
aerodynamic  heating  conditions  can  be  defined. 


V.  Test  Data  and  Analysis 

The  facility  in  which  was  gathered  the  majority 
of  the  data  to  be  reported  in  this  paper  was  the 
1.  5  MW  arc  heater  test  unit  of  the  Propulsion 
Wind  Tunnel  Facility  at  AEDC  in  Tullahoma,  Tenn. 
This  facility  uses  the  Linde  N-4000  heater  with  an 
operational  chamber  pressure  of  100  atmospheres 
maximum.  The  arc  heater  has  two  modes  of 
operation  in  terms  of  the  cathode.  The  short 
cathode,  26  inches  in  length,  gives  a  bulk  enthalpy 
of  3000  Btu  per  lb,  and  the  long  cathode,  34  inches 
in  length,  gives  2600  Btu  per  lb.  Three  nozzles 
were  used  during  this  test  program;  the  first  is  a 
Mach  2.  3  nozzle  and  yields  stagnation  pressure  on 
the  model  of  55  atmospheres  at  a  chamber  pres¬ 
sure  of  100  atmospheres.  The  second,  a  Mach 
2.0  nozzle,  yields  70  atmospheres  on  the  model, 
and  the  third,  a  Mach  1.  6  nozzle,  yields  85  atmos¬ 
pheres  on  the  model.  A  full  description  of  the 
facility  and  its  capability  is  given  in  Reference  42. 

Conditions 

The  hemisphere-cone  models  were  tested  in 
the  free  jet  configuration.  The  nozzle  expansion 
ratios  and  chamber  pressures  were  such  that  the 
free  jet  was  always  supersonic  and  highly  under¬ 
expanded.  The  characterization  of  the  test  environ¬ 
ment  involved  definition  of  the  model  flowfield, 
and  model  stagnation  enthalpy  and  pressure.  The 
jet  flow  field  calibration  is  described  in  Refer¬ 
ences  42  and  43,  and  thus  only  a  brief  summary 
will  be  given  here.  The  high  heating  rate  makes 
the  design  of  instrumentation  with  sufficient 
response  and  longevity  for  making  accurate  mea¬ 
surements  difficult.  For  this  reason,  measure¬ 
ment  of  heat  flux  and  of  other  phenomena  such  as 
turbulence,  pressure,  and  enthalpy  profiles  has 
not  been  attempted.  It  was  found  that  the  measured 


values  of  heat  flux  did  not  agree  with  those  predicted 
by  Fay-Riddell  (^4)  theory  and  any  of  the  above 
unmeasured  phenomena  could  contribute  to  this 
disagreement. 

It  is  characteristic  of  the  Linde  N-4000  arc 
heater  to  produce  a  flow  with  a  high  enthalpy  core; 
i.  e.  ,  the  centerline  stagnation  enthalpy  of  the  jet 
is  significantly  greater  than  the  average  bulk 
enthalpy  of  the  stream.  Since  insufficient  evalu¬ 
ation  of  the  turbulence  or  pressure  gradients  exists 
to  allow  quantitative  assessment,  the  approach  in 
this  work  has  been  to  assume  that  the  enthalpy 
profile  is  the  source  of  the  above-noted  discrepancy. 

To  determine  the  magnitude  of  the  centerline 
stagnation  enthalpy  of  the  jet,  heat  flux  measure¬ 
ments  were  made  using  high  response  transient 
calorimeters.  The  apparent  centerline  stagnation 
enthalpy  of  the  jet  was  then  calculated  using  the 
stagnation  point  heating  analysis,  of  Fay  and 
Riddell,  Figure  23  shows  the  correlation  of  appar¬ 
ent  centerline  stagnation  enthalpy  for  the  arc 
heater  configurations  used. 


Figun  23.  Ratio  of  Jat  Ointirtint  Stagnation  Enthaipy  to 
Bulk  Enthalpy 


Pertinent  arc  heater  parameters  were  mea¬ 
sured  during  each  run  to  allow  a  heat  balance 
determination  of  the  average  total  (bulk)  enthalpy 
of  the  flow.  The  total  pressure  of  the  flow  was 
measured  directly.  The  model  flowfield  could, 
therefore,  be  defined  for  each  arc  heater  firing 
period.  Ablation  model  data  included  surface 
brightness  temperature  measurements  and  high¬ 
speed  (400  to  600  frames /second)  motion  picture 
coverage.  Appropriate  lenses  and  filters  were 
used  on  the  cameras  to  provide  a  large  and  clear 
representation  of  the  ablating  surface.  All 
recorded  data  displays  were  correlated  by  an 
independent  timing  device.  The  teat  measure¬ 
ments  and  measurement  techniques  for  the  test 
program  reported  herein  are  given  in  Reference  45. 

Data  Reduction  and  Analysis 

A  total  of  140  graphite  models  were  tested  at 
AEDC,  including  fine-  and  coarse-grain  graphites. 
All  of  these  models  have  rather  small  nose  radii 
either  0.  18  or  0.25  in.  Surface  temperature  mea¬ 
surements  were  taken;  however,  difficulties  were 
encountered.  The  problem  lies  not  in  the  equip¬ 
ment  available,  but  in  the  position  in  which  this 
equipment  must  be  placed  in  order  to  see  the  sur¬ 
face  of  the  model.  The  model  proximity  to  the  jet 
exit  is  such  that  the  pyrometers  had  to  be  placed  at 
90  degrees  to  the  model  centerline.  Consequently, 
the  pyrometer  readings  reflected  an  average  tem¬ 
perature  over  the  periphery  of  the  model  hemis¬ 
phere  nose  and  shoulder  and  are  unacceptable  for 
the  stagnation  point  correlations  presented  here. 


The  recession  rates  have  been  determined 
from  the  motion  picture  coverage  in  order  to  obtain 
recession  measurements  after  model  warmup  and 
before  the  model  receded  beyond  the  point  where 
the  total  pressure  began  decreasing.  A  typical 
plot  of  data  is  shown  in  Figure  24.  Experimental 
recession  data  in  film  reader  units  versus  time 
from  model  injection  are  shown.  The  points  rep¬ 
resent  individual  measurements  made  from  frames 
of  the  movie  film.  After  a  short  warmup  period 
the  model  evidenced  linear  recession  and  was 
apparently  continuous.  This  is  important,  the 
linear  recession  indicated  shows  that  the  model 
rapidly  approached  steady  state  and  that  no  major 
observable  spallations  occurred.  The  model 
recession  rate  determined  is  the  slope  of  the  linear 
portion  of  the  data.  As  indicated,  recession  rates 
appear  to  decrease  as  the  model  recedes.  This  is 
caused  by  the  total  pressure  decrease  after  passing 
from  the  "cone  of  silence."  In  reducing  the  data, 
recession  rates  are  determined  from  the  linear 
portion  of  the  experimental  data,  thus  avoiding  the 
influence  of  nose  blunting  and  decay  in  heat  rate 
and  pressure  on  recession  rates  as  the  model 
recedes. 


Figun  24.  Typieil  Surfm  RiCMion  Hiitoiy  from  Film 
Rsiding 

A  typical  block  of  data  reduced  in  this  manner 
is  shown  in  Figure  25;  experimental  recession  data 
are  plotted  against  apparent  stagnation  enthalpy. 
Two  sets  of  data  are  given  as  typical.  Shown  are 
data  taken  in  the  stagnation  pressure  range  of 
83.  5  to  88.  6  atm  and  data  gathered  between  18.  6 
and  23.  1  atm.  Each  set  of  data  was  correlated 
using  a  least  squares  curve  fit  of  the  data  points. 
Results  similar  to  these  have  been  obtained  for 
nominal  stagnation  pressures  of  55  and  70  atm.  In 
addition,  recession  data  in  the  6.  5  to  9.8  atm  range 
have  been  determined  using  data  from  Reference  46. 
The  composite  of  the  recession  data  enables  con¬ 
struction  of  the  experimental  recession  rate  map 
shown  in  Figure  26.  Here  experimental  recession 
rate  is  shown  in  terms  of  stagnation  pressure  and 
stagnation  enthalpy.  Figure  26  presents  the  total 
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picture  of  the  dependency  of  the  experimental 
recession  rate  data  upon  the  environmental  condi¬ 
tions,  stagnation  enthalpy,  and  pressure.  The  data 
appear  to  be  linear  functions  of  pressure  and 
enthalpy. 


Fi|un  25.  Exparimsotsl  RscMion  Ran  Data 


Fifuia  26.  Exparimantil  Racaaakin  Rata  Data  Map 


Construction  of  the  experimental  recession 
data  map  enables  a  comparison  of  the  calculated 
thermochemical  recession  to  the  observed  experi¬ 
mental  recession  data.  Figure  27  presents  the 
ratio  of  experimental  and  thermochemical  mass 
fluxes  for  the  test  environmental  spectrum. 


Fi|ura  27.  Ratio  of  Exparimantal  and  Tfwrmochsmical  Maai 
Fhixai 


Thcrmochemical  recession  rates,  including  the 
uncertainty  band  in  the  carbon  vapor  properties, 
were  given  in  Figure  ll  for  this  environmental 
spectrum.  For  clarity,  this  uncertainty  in  ther- 
mochemical  recession  rates  as  affects  the  ratio 
given  in  Figure  27  is  omitted  (the  ratio  is  affected 
by  less  than  10%).  As  indicated  in  Figure  27  the 
deviation  from  thcrmochemical  recession  is 
marked,  being  affected  by  both  stagnation  pressure 
and  enthalpy.  The  experimental  mass  flux  is  about 


a  factor  of  2  greater  than  that  predicted  by  thermo¬ 
chemical  erosion  alone  at  100  atm  pressure.  Com¬ 
plementary  results  are  given  in  Figure  28  depicting 
the  predicted  surface  temperature  as  affected  by 
loss  of  solid  from  the  heated  surface.  Noteworthy, 
is  the  rather  small  decrease  in  surface  tempera¬ 
ture  for  even  large  fractions  of  solid  eroded.  Sur¬ 
face  temperature  is  most  affected  by  solid  loss  at 
the  low  stagnation  enthalpy;  the  high  stagnation 
enthalpy  somewhat  masks  the  less  efficient 
removal  of  material  in  the  solid  state. 


Fi|ura28.  InfkMKS  of  Solid  Friction  Eroded  on  Surfaco 
Timpontun 


Figure  29  shows  the  departure  from  thermo¬ 
chemical  erosion  in  more  meaningful  terms.  From 
the  composite  of  experimental  data  (Figure  26)  and 
numerical  results,  including  solid  removal,  the 
fraction  of  solid  eroded  has  been  calculated  and  is 
depicted  in  Figure  29.  Coincident  with  the  results 
given  in  Figure  27  the  threshold  stagnation 
pressure  at  which  loss  of  particulate  solid  initiates 
is  dependent  on  stagnation  enthalpy  (reflecting  con¬ 
vective  heat  rate  at  constant  pressure).  The 
results  shown  are  consistent  with  those  expected, 
increased  stagnation  enthalpy  results  in  higher 
surface  temperatures;  hence,  greater  subsurface 
loss  of  binder  material  reducing  binding  of  filler 
particles  within  the  graphite  matrix.  Increasing 
pressure  above  the  threshold  value  causes  a  large 
increase  in  the  solid  fraction  eroded  over  a  fairly 
narrow  pressure  range.  The  Increment  in  solid 
eroded  thereafter  decreases  as  the  full  pressure 
range  is  covered. 


Figun  29.  Fraction  of  SoNd  Erodid  from  Expirbnontal  Dott 
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To  interpret  the  results  in  Figure  29  in  terms 
of  the  contents  of  this  paper  reference  is  made  to 
the  postulated  preferential  erosion  mechanism  pre¬ 
viously  discussed.  First,  it  was  observed 
(Figure  28)  that  the  effect  of  pressure  on  surface 
temperature  is  small  at  the  respective  values  of 
stagnation  enthalpy  (within  the  pressure  range 
covered).  Thus,  each  value  of  stagnation  enthalpy 
can  be  viewed  in  terms  of  nearly  constant  surface 
temperature  for  the  range  of  solid  fraction  eroded 
evidenced  in  Figure  29.  Consider  then  the  role  of 
stagnation  pressure,  the  normal  force  acting  on  the 
heated  surface.  As  the  surface  recedes,  filler 
particles  of  various  siee  are  continuously  exposed. 
Loss  of  binder  below  the  heated  surface  and  by 
direct  thermochemical  erosion  at  the  surface  tends 
to  free  the  otherwise  constrained  and  somewhat 
interlocked  particles.  The  larger  filler  particles 
are  then  rather  easily  blown  free  by  even  moderate 
pressures,  commencing  with  the  higher  tempera¬ 
ture  exposure  conditions  (i.  e. ,  stagnation  enthalpy) 
with  increased  pressure  acting  to  free  the  smaller 
particles.  Also,  minute  amounts  of  binder  leave 
the  surface  attached  to  the  dislodged  filler  particles.. 
As  indicated  in  Figure  29  the  effect  of  enthalpy 
(hence,  surface  temperature)  at  the  higher  pres¬ 
sures  appears  to  diminish  and  pressure  dominates 
the  loss  of  particulate  solid.  Interestingly,  the 
fraction  of  solid  eroded  appears  to  be  approaching 
an  asymptotic  value  in  the  range  0,  6  to  0.  8. 
Significantly,  the  nominal  mass  fraction  of  filler 
particles  (about  0.  75)  falls  within  this  range. 

Finally,  the  effect  of  aerodynamic  shear  and  pres¬ 
sure  gradients  typical  of  positions  removed  from 
the  stagnation  point  is  to  induce  increased  particle 
loss.  Thus,  one  is  faced  with  loss  of  particulate 
graphite  over  a  rather  broad  range  of  aerodynamic 
conditions,  perhaps  even  extending  into  the  rela¬ 
tively  moderate  super-orbital  entry  environment. 


VI.  Conclusions 

A  model  incorporating  the  allotropic  features 
of  bulk  graphite  has  been  postulated  and  used  to 
describe  the  continuous  removal  of  particulate 
graphite  from  the  heated  surface.  Although  hypo¬ 
thetical  in  its  formation,  sufficient  laboratory  data 
and  pre-  and  post-test  arc-jet  data  are  advanced 
to  establish  the  creditability  of  removal  of  filler 
particles  in  the  solid  state.  The  key  role  of  the 
different  forms  of  carbon  in  graphite  on  surface 
erosion  has  been  established.  Thermodynamic 
properties  of  the  linear  and  cyclic  molecules  C3  to 
C30  were  calculated  by  quantum  and  statistical 
mechanical  techniques.  The  resonance  energy  of 
each  molecule  was  computed  from  the  respective 
electronic  wave  functions,  and  the  corresponding 
vibrational  frequencies  were  determined  by  a 
normal  coordinate  analysis  of  the  molecular  struc¬ 
ture.  The  theoretical  procedures  are  self- 
consistent,  resort  to  the  empirical  increments 
adopted  by  Pitser  and  dementi  and  used  by  others 
has  been  avoided.  The  agreement  of  the  calculated 
vapor  pressure  with  the  triple-point  and  sublima¬ 
tion  data  is  excellent.  By  purely  analytical  means, 
the  uncertainty  in  the  free  energy  of  formation  has 
been  established;  the  effect  on  calculated  thermo¬ 
chemical  surface  recession  was  found  to  be  small. 
Significantly,  this  result  has  enabled  identification 
of  the  solid  fraction  of  graphite  eroded  in  high- 
pressure  ground  tests  with  confidence. 


To  establish  the  solid  fraction  of  particulate 
graphite  eroded  ablation  data  encompassing  fine  to 
coarse  grain  graphites  have  been  correlated.  The 
loss  of  material  in  the  solid  state  has  been  deter¬ 
mined  over  a  range'  of  environmental  conditions  (in 
terms  of  stagnation  pressure  and  enthalpy)  evi¬ 
dencing  large  fractions  of  solid  removal.  Loss  of 
particulate  solid  initiates  at  rather  low  stagnation 
pressures;  the  actual  threshold  is  dependent  on 
heating  rate  as  well  as  pressure.  The  mechanism 
by  which  loss  of  particulate  graphite  occurs 
appears  to  be  by  preferential  surface  recession  and 
subsurface  loss  of  binder  phase,  freeing  the  filler 
particles  from  the  binder  matrix.  Filler  particles 
are  then  easily  dislodged  from  the  surface  by  even 
moderate  forces  (i.  e. ,  pressure,  shear).  Surpris¬ 
ingly,  the  loss  of  solid  at  the  heated  surface  does 
not  appreciably  depress  surface  temperatures, 
even  for  rather  large  fractions  of  solid  loss. 

The  results  presented  in  this  paper  open 
several  interesting  avenues  for  pursuit.  Clearly, 
the  role  of  graphite  microstructure  in  dictating 
surface  erosion  cannot  be  dismissed.  Closely 
related  to  the  microstructure  of  the  binder  and 
filler  phases  are  the  thermophysical,  mechanical, 
and  chemical  properties  which  individually,  or  in 
combination,  affect  ablation  performance.  Future 
analyses  employing  statistical  distributions  of 
particle  size,  void  size  and  spacing,  and  account¬ 
ing  for  the  temperature-time  course  of  the  crystal¬ 
lite  growth  will  be  necessary  before  a  completely 
analytical  performance  description  of  graphite  can 
be  achieved.  Finally,  new  graphites  overcoming 
the  evident  dissimilar  characteristics  of  the  con¬ 
stituent  materials  appear  to  afford  a  means  of 
reducing  or  alleviating  preferential  erosion,  hence, 
loss  of  solid  material. 
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